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When a particle is driven with an ac force over an asymmetric potential, it can undergo a ratchet 
effect that produces a net dc motion of the particle. Ratchet effects have been observed in numerous 
systems such as superconducting vortices on asymmetric pinning substrates. Skyrmions, stable 
topological spin textures with particle-like properties, have many similarities to vortices but their 
behavior is strongly influenced by non-disspative effects arising from a Magnus term in their equation 
of motion. We show using numerical simulations that pronounced ratchet effects can occur for ac 
driven skyrmions moving over asymmetric quasi-one-dimensional substrates. We find a new type 
of ratchet effect called a Magnus-induced transverse ratchet that arises when the ac driving force 
is applied perpendicular rather than parallel to the asymmetry direction of the substrate. This 
transverse ratchet effect only occurs when the Magnus term is finite, and the threshold ac amplitude 
needed to induce it decreases as the Magnus term becomes more prominent. Ratcheting skyrmions 
follow ordered orbits in which the net displacement parallel to the substrate asymmetry direction 
is quantized. Skyrmion ratchets represent a new ac current-based method for controlling skyrmion 
positions and motion for spintronic applications. 


Ratchet effects can arise when particles placed in 
an asymmetric environment are subjected to suitable 
nonequilbrium conditions such as an externally applied 
ac drive, which produces a net dc motion of the particles^. 
Ratchet effects have been realized for colloidal systems^, 
cold atoms in optical traps^, granular media on sawtooth 
substrates^, and the motion of cells crawling on patterned 
asymmetric substrates^, and they can be exploited to 
create devices such as shift registers for domain walls 
moving in an asymmetric substrate^. There has been in¬ 
tense study of ratchet effects for magnetic flux vortices 
in type-II superconductors interacting with quasi-one- 
dimensional (qlD) or two-dimensional (2D) asymmetric 
pinning substrates. The vortices can be effectively de¬ 
scribed as overdamped particles moving over an asym¬ 
metric substrate, and under ratcheting conditions the 
application of an ac driving current produces a net dc 
flux flow£~— . The simplest vortex ratchet geometry was 
first proposed by Lee et al. 7 , where an effectively 2D as¬ 
sembly of vortices was driven over a qlD asymmetrically 
modulated substrate using an ac driving force oriented 
along the direction of the substrate asymmetry. 

Skyrmions in chiral magnets, initially discovered 
in MnSii^i and subsequently identified in numerous 
other materials at lowiA— and room temperatures^ - —, 
have many similarities to superconducting vortices. 
Skyrmions are spin textures forming topologically sta¬ 
ble particle-like objects^ that can be set into motion 
by the application of a spin-polarized currenl— - — . As 
a function of the external driving, skyrmions can ex¬ 
hibit a depinning transition similar to that found for 
vortices, and from transport measurements it is possi¬ 
ble to construct skyrmion velocity versus applied force 
curve s 22 ' 24- —. It has been shown that the dynamics of 
skyrmions interacting with pinning can be captured by 
an effective particle equation of motion, which produces 


pinning-depinning and transport properties that agree 
with those obtained using a continuum-based model^, 
just as a particle-based description for vortices can ef¬ 
fectively capture the vortex dynamics in the presence of 
pinning sites— Due to their size scale and the low cur¬ 
rents required to move them, skyrmions show tremen¬ 
dous promise for applications in spintronics such as race¬ 
track type memory device s 30 ' 31 originally proposed for 
domain walls^, and it may be possible to create skyrmion 
logic devices similar to those that harness magnetic do¬ 
main walls— In order to realize such skyrmion applica¬ 
tions, new methods must be developed to precisely con¬ 
trol skyrmion positions and motion. 

The most straightforward approach to a skyrmion 
ratchet is to utilize the same types of asymmetric 
substrates known to produce ratchet effects in vor¬ 
tex systems. There are, however, important differ¬ 
ences between skyrmion and vortex dynamics due to the 
strong non-dissipative Magnus component of skyrmion 
motion^ ' — il 22 ' 24 ' 2 -' ' 2 ^ that is absent in the vortex sys¬ 
tem. The Magnus term produces a skyrmion velocity 
contribution that is perpendicular to the direction of an 
applied external force. In the presence of pinning, the 
Magnus term reduces the effectiveness of the pinning by 
causing the skyrmions to deflect around the edges of at¬ 
tractive pinning site s—'— '— rather than being captured 
by the pinning sites as overdamped vortices are. The 
Magnus term also produces complex winding orbits for 
skyrmions moving in confined regions^ or through pin¬ 
ning site s 24 ' 26 ' 27 ' 35- — . The ratio of the strength of the 
Magnus term a m to the dissipative term ay can be 10 
or higher for skyrmion systems— In contrast, although 
it is possible for a Magnus effect to appear for vortices 
in superconductors, it is generally very weak so that the 
vortex dynamics is dominated by dissipation—. A key 
question is how the Magnus force could impact possible 
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FIG. 1: Asymmetric substrate geometry: A schematic 
of the sample geometry where a skyrmion is placed on a 
qlD asymmetric periodic substrate. An ac driving force 
F ac cos (uit) is applied parallel (_Fjf c ) or perpendicular ( F “ c ) 
to the substrate asymmetry direction. 
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ratchet effects for skyrmions, and whether new types of 
ratchet phenomena can be realized for skyrmion systems 
that are not accessible in overdamped systems^. 

In this work we numerically examine a 2D assembly of 
skyrmions interacting with a qlD asymmetric periodic 
substrate shown schematically in Fig. 1. An ac drive is 
applied either parallel (Fff c ) or perpendicular (F“ c ) to 
the substrate asymmetry direction, which runs along the 
x axis. In the overdamped limit, no ratchet effect ap¬ 
pears for F“ c ; however, we find that the Magnus effect 
produces a novel ratchet effect for F“ c when it curves 
the skyrmion trajectories into the asymmetry or x di¬ 
rection. We model the skyrmion dynamics in a sample 
with periodic boundary conditions using a particle based 
descriptio n 26 : 37 and vary the ratio of the Magnus to dis¬ 
sipative dynamic terms, the amplitude and frequency of 
the ac drive, and the strength of the substrate. 

Simulation 

The equation of motion for a single skyrmion with ve¬ 
locity v moving in the x — y plane is 

«dv + a m z x v = F sp + Fjf^ (1) 

The damping term ad keeps the skyrmion velocity 
aligned with the direction of the net external force, while 
the Magnus term a m rotates the velocity toward the di¬ 
rection perpendicular to the net external forces. To ex¬ 
amine the role of the Magnus term we vary the relative 
strength a m /ad of the Magnus term to the dissipative 
term under the constraint a^ + o? m = 1, which maintains 
a constant magnitude of the skyrmion velocity. The sub¬ 
strate force F sp = —V17(a;)x arises from a ratchet po¬ 
tential 

U(x) = f7o[sin(27nc/a) + 0.25 sin(47rx/a)], (2) 

where a is the periodicity of the substrate and we de¬ 
fine the strength of the substrate to be A p = 2nUo/a. 


FIG. 2: Parallel and perpendicular ratchet veloci¬ 
ties for parallel drive: The skyrmion velocity parallel (F)|| 
(red) and perpendicular (F)j_ (blue) to the substrate asym¬ 
metry vs F{/ c for a substrate strength of A p = 1.5. (a) At 
Om/ccd = 0.0, the overdamped limit, only (F)|| is nonzero 
and shows quantized steps with (F)ii = n, where n is an inte¬ 
ger. (b) At am/ad, = 1.6, there is quantized motion in both 
the parallel and perpendicular directions, (c) am/ad. = 4.58. 
(d) am/ad = 7.85. The points labeled B and C indicate where 
the skyrmion orbits in Fig. 3(b,c) were obtained. 


The ac driving term is either Fjj c = Fj“ c cos(wt)x or 
F“ c = F“ c cos(wt)y. We measure the time-averaged 
skyrmion velocities parallel (F)|| = 27 t(v • x)/wa or per¬ 
pendicular (F)j_ = 27 t(v • y )/wa to the substrate asym¬ 
metry direction as we vary the ac amplitude, substrate 
strength, or u>. Here we focus primarily on samples 
with A p = 1.5 and oj = 2.5 x 10~ 6 inverse simulation 
time steps, and we consider the sparse limit of a sin¬ 
gle skyrmion which can also apply to regimes in which 
skyrmion-skyrmion interactions are negligible. 


Results and Discussion 

A. ac Driving Parallel to Substrate Asymmetry 

We first apply the ac driving force parallel to the asym¬ 
metry direction, a configuration previously studied for 
the same type of substrate in vortex systemsid^,. In 
Fig. 2 we plot <F)|| and (F)j_ versus F^ c for samples 
with A p = 1.5, so that the maximum magnitude of the 
substrate force is F^ ax = 2.0 in the negative x direction 
and F^ ax = 1.0 in the positive x direction. The over¬ 
damped limit a m /ad = 0.0 appears in Fig. 2(a), where 
{V)± = 0.0 and there is a ratchet effect only along the 
parallel or cc-direction for F^ c > 1.25. The skyrmion 
velocity is quantized and forms a series of steps with 
(F)|| = n where n is an integer. On the highest step 
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in Fig. 2(a) near F^ c = 2.3, n = 7. The quantization in¬ 
dicates that during one ac driving period, the skyrmion 
moves a net distance of na in the parallel direction. 
For very large values of F^ c , the velocity quantization 
is lost and the parallel ratchet effect gradually dimin¬ 
ishes. In Fig. 2(b) we plot (V)|| and (V)_l for a sample 
with a m /ad = 1.6, showing a nonzero skyrmion veloc¬ 
ity component for both the parallel and perpendicular 
directions. There is again a quantization of the parallel 
velocity {V) ||, and the Magnus term transfers this quanti¬ 
zation into the perpendicular direction even though there 
is no periodicity of the substrate along the y direction. 
Thus, we find {V)x = Tia m /ad■ For F^ c > 5.0, we 
find windows in which the ratcheting effect is lost and 
(Id)|| = (y)± = 0. As a m /ad increases, the maximum 
value of n decreases, so that at a m /ad = 4.58 in Fig. 2(c), 
only the n = 1 step appears. At the same time, the 
width of the non-ratcheting windows increases, as shown 
in Fig. 2(d) for a m /ad = 7.858. These results show that 
while a ratchet effect does occur for skyrmion systems, 
an increase in the Magnus term decreases the range of ac 
drives over which the ratchet effect can be observed. 

By lowering the ac frequency ui or increasing the sub¬ 
strate strength A p . it is possible to increase the range and 
magnitude of the skyrmion ratchet effect. In Fig. 3(b,c) 
we show representative skyrmion orbits for the system in 
Fig. 2(b,c). Figure 3(b) shows the orbits at a m /ad = 1.6 
for F|“ c = 2.2, corresponding to the n = 3 step in 
Fig. 2(b). Here the skyrmion moves along straight lines 
oriented at an angle to the direction of the applied ac 
drive. During each ac cycle, the skyrmion first translates 
a distance 3 a in the positive x direction and then travels 
in the reverse direction for a distance of a/2 before strik¬ 
ing the top of the potential barrier, which it cannot over¬ 
come in the reverse x direction for this magnitude of ac 
drive. It repeats this motion in each cycle. In Fig. 3(c) we 
show the skyrmion orbit at a m /ad = 4.58 for F’|“ c = 2.0, 
corresponding to the n = 1 orbit in Fig. 2(c). The angle 
the skyrmion motion makes with the ac driving direc¬ 
tion is much steeper, and the skyrmion is displaced a net 
distance of a in the x direction during each ac driving 
cycle. 

In order to clarify the evolution of the ratchet phases 
as a function of a m /ad and F^ c , in Fig. 3(a) we plot the 
ratchet envelopes for the n = 1 to n = 6 steps. For the 
n = 1 case we show, in blue, all of the regions in the 
F|“ c — a m /ad plane where n = 1 steps appear. For n = 2 
to n = 6, to prevent overcrowding of the graph, instead 
of plotting all ratchet regions we show only the inner step 
envelope for each n = rii, obtained from simulations in 
which Fj“ c is swept up from zero, and defined to extend 
from the drive at which n first reaches rq to the drive at 
which n first drops below rq . We find that the minimum 
F|“ c required to induce a ratchet effect increases linearly 
with increasing a m /ctd , and that the ratcheting regions 
form a series of tongue features (shown for n = 1 only). 
The ratchet effects persist up to and beyond a m /ad = 10. 



FIG. 3: Ratchet region envelopes and example ratchet 
orbits for parallel ac driving: (a) A plot of F{/ c vs a m /ctd 
indicating all n = 1 (blue) ratcheting regions, and the inner 
envelopes of the n = 2 (orange), 3 (purple), 4 (green), 5 
(red) and 6 (yellow) ratcheting regions. The boundaries of the 
inner envelopes are defined as the point at which the ratchet 
velocity first reaches the value n followed by the point at 
which it first drops below n. The y-axis marked A and the 
dashed lines marked B — D are the values of am/ad for which 
the velocity curves in Fig. 2 were obtained. The minimum ac 
amplitude required to produce a ratchet effect increases with 
increasing am/ad■ (b,c) Skyrmion trajectory images. White 
(green) regions are high (low) areas of the substrate potential. 
Thin black lines indicate forward motion of the skyrmion and 
thick purple lines indicate backward motion. The red circle 
is the skyrmion. (b) The n = 3 orbit at am/ad = 1.6 for 
F|“ c = 2.2 from Fig. 2(b). The skyrmion moves in straight 
lines and translates 3 a in the x direction during one ac drive 
cycle, (c) The n = 1 orbit at a m /ad = 4.58 for Fj“ c = 2.0 
from Fig. 2(c), where the skyrmion translates a distance a in 
the x direction during every ac drive cycle. 


For a m /ad = 0, the strongest ratchet effects with the 
largest values of n occur near Fj“ c = 2.3, corresponding 
to a driving force that is slightly higher than the max¬ 
imum force exerted on the skyrmion by the substrate 
when it is moving with a negative x velocity component. 
As F|“ c increases above this value, the skyrmion can slip 
backward over more than one substrate plaquette during 
each ac cycle, limiting its net forward progress. 
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FIG. 4: Parallel and perpendicular ratchet veloci¬ 
ties for perpendicular drive: The skyrmion velocity par¬ 
allel (V’)n (red) and perpendicular (P)x (blue) to the sub¬ 
strate asymmetry vs F/f for samples with A p = 1.5. (a) At 
am/ctd = 0.855, ratcheting occurs in both the parallel and 
perpendicular directions for E“ c > 1.5, and (F)n = n with 
integer n. Inset: At a m /ad = 0, there is no ratchet effect, 
(b) a m /ad = 4.0. (c) am/ad = 7.018. Inset: A blowup of 
(V’)ll from the main panel at the third n = 1 step showing 
that additional fractional steps appear at velocity values such 
as n/m = 1/2. 


B. ac Driving Perpendicular to Substrate 
Asymmetry 



FIG. 5: Example ratchet orbits for perpendicular 

ac driving: Skyrmion trajectory images from the system 
in Fig. 4. White (green) regions are high (low) areas of the 
substrate potential. Black lines indicate the motion of the 
skyrmion; the red circle is the skyrmion. (a) The n = 2 step 
at a m /ad = 0.855 for E/ c = 2.0. (b) The n = 1 step at 
Om/od = 4.0 for F“ c = 0.7. (c) A non-ratcheting orbit at 
a m /ad = 4.0 for E“ c = 0.76. (d) For higher ac drives the or¬ 
bits become more complicated, as shown here for a m /ad = 4.0 
at Fl c = 1.97. 


In Fig. 4 we plot (P)|| and (V)j_ versus F“ c for A p = 
1.5. The inset in Fig. 4(a) shows that for a m /ad = 0, 
(V')|| = (V)j. = 0, indicating that there is no ratchet ef¬ 
fect in either direction. When the Magnus term is finite, 
as in Fig. 4(a) where a m /ad = 0.855, pronounced ratch¬ 
eting occurs in both the parallel and perpendicular direc¬ 
tions. We describe this as a Magnus-induced transverse 
ratchet effect. Here (V)|| is quantized at integer values, 
with the largest step in Fig. 4(a) falling at n = 5, and 
{V)±/{V)\\ = am/ad . As a m /a d increases, Fig. 4(b,c) 
shows that the maximum value of F“ c at which ratchet¬ 
ing occurs decreases, as does the maximum value of n and 
the widths of the ratcheting windows. In addition to the 
integer steps in (V)||, we also observe fractional steps, as 
shown in the inset of Fig. 4(c) for a m /ad = 7.018, where 
we highlight the third n = 1 step. On the upper n = 1 
step edge there is a plateau at (V)|| =1/2 and additional 
steps at (V’)n = 3/4, 2/3, 1/3, and 1/4, which are accom¬ 
panied by even smaller steps that form a devil’s staircase 
structure^. For higher values of F// c , smaller plateaus 
appear at rational fractional velocity values (V')|| = n/m 
with integer n and m, while full locking to the integer 
n = 1 step no longer occurs. The n/m rational lock¬ 
ing steps also occur for parallel ac driving but are much 
weaker than for perpendicular ac driving. 


The skyrmion orbits for perpendicular ac driving differ 
markedly from those that appear under parallel ac driv¬ 
ing, since the Magnus term induces a velocity compo¬ 
nent perpendicular to the direction of the external driv¬ 
ing force. As a result, perpendicular ac drives induce 
parallel skyrmion motion that interacts with the sub¬ 
strate asymmetry. Figure 5(a) shows a skyrmion orbit 
on the (V)|| = 2 step for F“ c = 2.0 and a m /ad = 0.855 
from Fig. 4(a). The skyrmion no longer moves in straight 
lines, as it did for the parallel ac drive in Fig. 3. Instead 
it follows a complicated 2D trajectory, translating a net 
distance of 2a in the x direction during every ac drive 
cycle. In Fig. 5(b) we plot an n = 1 orbit at F“ c = 0.7 
for a m /ad = 4.0 from Fig. 4(b), where there is a sim¬ 
ilar orbit structure in which the skyrmion translates a 
net distance of a in the x direction in each ac drive cycle. 
Figure 5(c) shows a non-ratcheting orbit for a m /ctd = 4.0 
at F r / C = 0.75, where the skyrmion traces out a compli¬ 
cated periodic closed cycle. For higher ac amplitudes, the 
ratcheting orbits become increasingly intricate, as shown 
in Fig. 5(d) for a m /ud = 4.0 at F“ c = 1.97. 

In Fig. 6(a) we show the full n = 1 ratchet regions and 
the inner envelopes of the n = 2 to 5 ratchet regions as 
a function of F“ c and a m /ad■ For low values of a m /ad-, 
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FIG. 6: Ratchet region envelopes for perpendicular 
ac driving: (a) A plot of F“ c vs am/ad indicating all n = 1 
(blue) ratcheting regions, and the inner envelopes of the n — 2 
(orange), 3 (purple), 4 (green), and 5 (red) ratcheting regions. 
The dashed lines marked A — C are the values of o m /ad for 
which the velocity curves in Fig. 4 were obtained. There is no 
ratchet effect at the overdamped limit of otm/ud = 0, and the 
minimum ac amplitude required to produce a ratchet effect 
decreases with increasing a m /ad ■ (b) A plot of F“ c vs A p for 
otm/otd = 4.92 indicating the ratcheting regions with n = 1 
(blue), 2 (orange), 3 (purple), and 4 (green). As the substrate 
strength increases, stronger ratchet effects appear. 


there is no ratcheting since the Magnus term is too weak 
to bend the skyrmion trajectories far enough to run along 
the direction parallel to the asymmetry. The minimum ac 
amplitude required to generate a ratchet effect decreases 
as the strength of the Magnus term increases, opposite to 
what we found for parallel ac driving in Fig. 3(a). Here, 
the skyrmion trajectories are more strongly curved into 
the parallel direction as a m /ctd increases, so that a lower 
ac amplitude is needed to push the skyrmions over the 
substrate potential barriers. In contrast, for a parallel ac 
drive the Magnus term curves the skyrmion trajectories 
out of the parallel direction, so that a larger ac ampli¬ 
tude must be applied for the skyrmions to hop over the 
substrate maxima. We also find that as a m /ad increases, 
the minimum ac force needed to produce a ratchet effect 
drops below the pinning force exerted by the substrate. 
This cannot occur in overdamped systems, and is an in¬ 
dication that the Magnus term induces some inertia-like 
behavior. The threshold ac force value F th for ratcheting 
to occur can be fit to F^' oc ( a m /ad)~ l for perpendicu¬ 
lar driving, while for parallel driving, Ft h oc a m /ad■ As 
otm/ocd increases, the extent of the regions where ratch¬ 
eting occurs decreases; however, by increasing the sub¬ 
strate strength, the regions of ratcheting broaden in ex¬ 
tent and the magnitude of the ratchet effect increases. 
In Fig. 6(b) we show the evolution of the n = 1, 2, 3 
and 4 ratchet regions as a function of F“ c and A p for 
fixed a m /ad = 4.92. Here, as A p increases, the mag¬ 
nitude of the ratchet effect increases, forming a series 
of tongues. The minimum ac force needed to induce a 
ratchet effect increases linearly with increasing A p . These 
results indicate that skyrmions can exhibit a new type 



FIG. 7: Parallel ratchet velocity heat maps for par¬ 
allel and perpendicular ac driving: The magnitude of 
(F)ll in samples with A p = 1.5 and am/ctd = 3.042. (a) As a 
function of F“ c and u> in units of 10~ 5 inverse simulation time 
steps, (F)|| has clear tongue structures, and its magnitude is 
non-monotonic for varied frequency, (b) (V) y as a function 
of F^ c and ta has similar features. 


of Magnus-induced ratchet effect that appears when the 
driving force is applied perpendicular to the asymmetry 
direction of the substrate. 


C. Frequency Dependence 

In Fig. 7(a) we plot a heat map of the parallel ratchet 
velocity (F)m for perpendicular ac driving as a function of 
F“ c and u in a sample with A p = 1.5 and a m /ad = 3.042. 
The ratchet effect extends over a wider range of F“ c at 
low ac driving frequencies, while for higher w the ratchet¬ 
ing regions form a series of tongues. The ratchet magni¬ 
tude is non-monotonic as a function of w, and in several 
regions it increases in magnitude with increasing ui, such 
as on the first tongue where the maximum ratcheting ef¬ 
fect effect occurs near w = 10~ 3 inverse simulation time 
steps before decreasing again as w increases further. We 
find very similar ratchet behaviors for parallel ac driving, 
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as shown in Fig. 7(b) where we plot (V)|| as a function 
of F|“ c and u and observe a series of tongues. 

Our results show that skyrmions are an ideal sys¬ 
tem in which to realize ratchet effects in the presence 
of asymmetric substrates. As expected, they undergo 
ratcheting when an ac drive is applied parallel to the 
substrate asymmetry direction; however, the skyrmions 
also exhibit a unique ratchet feature not found in over¬ 
damped systems, which is the transverse ratchet effect. 
Here, when the ac drive is applied perpendicular to the 
substrate asymmetry direction, the Magnus term curves 
the skyrmion orbits and drives them partially paral¬ 
lel to the asymmetry, generating ratchet motion. The 
asymmetric substrates we consider could be created us¬ 
ing methods similar to those employed to create quasi- 
one-dimensional asymmetric potentials in superconduc¬ 
tors and related systems, such as nano fabricated asym¬ 
metric thickness modulation, asymmetric regions of ra¬ 
diation damage, asymmetric doping, or blind holes ar¬ 
ranged in patterns containing density gradients. Since 


the skyrmion ratchet effects persist down to low frequen¬ 
cies, it should be possible to directly image the ratchet¬ 
ing motion, while for higher frequencies the existence of 
ratchet transport can be deduced from transport studies. 
Such ratchet effects offer a new method for controlling 
skyrmion motion that could be harnessed in skyrmion ap¬ 
plications. We also expect that skyrmions could exhibit 
a rich variety of other ratchet behaviors under different 
conditions, such as more complicated substrate geome¬ 
tries, use of asymmetric ac driving, or collective effects 
due to skyrmion-skyrmion interactions. 
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